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q. INTRODUCTION 


The ccmplexity involved in designing a control system 
cf reasonably low crder for dynamic systems cf substantially 
higher order is a subject about which much has been written. 
The list cf references at the end of this study include only 
a small ferticn of the material available concerning this 
subject. A new development has been intrcduced by Maurice 


F. Hutton and Bernard Friedland, [Ref 1], which provides a 


systematic and analytical approach eo obtaining 
reduced-crder approximations from higher-crder transfer 
fmurcticns. This method 1s referred to as the "Routh 


Approximaticr Methcd," and is based on an expansion that 
uses the Rceuth table of the original higher-order transfer 
function. Tke desire for low ocrder models to  Sirulate 
higher-order practical systems, such as electrical focwer 
plants, chemical processes, aircraft model designs and 
electrcnic cirfcwitry is obvious when dealing with such 
ccmplex systezs. The "Routh Approximation Method", is thus, 
an extremely powerful tool to utilize in dealing with these 


types of proklems. 


In ccmpariscn to other available methods presently 
used, the Kouth methcd has many advantages. The widely used 
“Dcminant Pcle'', approximation method , which is based on 
approximating a system by ted een Gg the poles, Or 
characteristic rocts nearest the imaginary axis, has the 
disadvantage that the roots of the characteristic equation 
must be found. For a very~high-order system, this is not a 
trivial task. The "Pad‘'e" approximation is tased on setting 
the numerator and denominator crders tc a desired value, the 
ccefficierts are then chosen so that the Taylor. series 


expansions cf the approximant and the original transfer 


EO 





function agreé€é in aS Many terms aS possible. This method 
preduces accurate results, however, it is limited in 
applicaticn tc single-input, single-output system analysis. 
In additicn, an unstable approximant may be obtained frem a 
stacle system since the approximant's poles depend on toth 


tke original eguation's numerator and denominator. 


The "Recuth Approximation Method," preserves stability 
if the criginal transfer function is stable. It provides an 
efficient reans of obtaining lower-order approximants for 
Multiple inputs or ottputs and is very adaptable to computer 


programming. 


This thesis is concerned with the development of a 
computer fregram which utilizes the "Routh Approxigation 
Method," to obtain lower-crder transfer functicns from 
higher-crder functions and compares the resulting éguations 
to the cridgiral eguations by graphically displaying the 
response cf each tc various input excitations. The prcegran 
was designed to provide output data which is useful in 
determining which degree of approximant is best suited to 
Simulate tke original higher-order equation. The crders 
available are the first through the fourth order reduced 
transfer functicns. The program is capable of recéucing 
transfer functions up to, and including the tenth orcer and, 
Without loss of generality, may be extended to handle an 


almest unlimited order. 


For iliustrative purposes, a seventh crder system was 
used, aS an exampie, to demonstrate the computer program's 
Capabilities. This particular system was chosen merely to 
demcnstrate the simplicity involved in utilizing the pregram 


and has ne physical relationship to any real systen. 


el: 


Peeve cone CE THE PROBL EM 


In e@ppreximating higher-crder systems by lower-order 
models, a linear system is desirable since the linear 
characteristic equations are less complicated than the 
non-linear eégquaticns. Therefore, this study is primarily 
ccncerned with linear, time invariant Characteristic 
equations, since the objective in develcfing the computer 
prcgram sas te illustrate a simple analytical approach to 
Obtaining lcwer crder characteristic eguations frem higher 


order systers. 


ite is obvicus that a reduced-order model cannot 
characterize a given system as accurately as a higher-order 
model. The validity of the lower-order model is based upon 
its degree of succeSS in approximating the higher order 
system in representing the characteristics of primary 


mimcerest. 


Interpreting the soluticn of a higher-order systen 
cften results in computational difficulties which are 
recuced Ey appropriate selection of a reduced-order 


apPEICcXimaticn. 


Ideally, a reduced-order model would approximate the 
higher-orcer system in both low and high frequency ranges. 
In dceing sc, scme accuracy is lost in compensating fcr the 
different responses of the systen eo Variations a 


frequency. 
The lew frequency model more closely approximates the 


higher-order system than a model composed cf both lew and 


high crder frequency characteristics. The procedure 


eZ 





utilized in this study places emphasis on the low freguency 


model. 


For unstable systems, the program described herin is 
still valid, however, the original higher-order systen must 
first be mecdified by shifting the imaginary axis prior to 
the apprcximation. Thus fOr an unstable transfer 
function,#(s), the equation must be changed tc H(sta), where 
a is chcser sufficiently large so that H (s) is 
asymptotically stable. This procedure is described in 
G@etadil itr f[kef t], pp 332. 


Le 





MeiereoneUrha PROGRAM CRITERIA 


A. GHNERAL 


The frogram, referred to as ROUTH1, waS written with 


the £collcwing criteria: 


ROUTH1 consists of less than 150K of storage and taxes 
less than 2 minutes of computer time. Pies. SC Gaon 
conserves efficiency, but also provides the user with the 


desired data at a minimum cost. 


2: Ease ¢ 


cel 


Ib 


use 
oa 


Input data reguired consists of eleven data cards for 
Maximum utilizaticn of the program. Emphasis is placed on 
the €ase cf using the program to obtain the desired results 


With minimum time expended on computer programing. 


3. Usefullness of output 


coy oad == 


a. Riometricst second, third and fourth order 
approximants to the original equation are printed in 
transfer function format. Both numerator and denominator 


coefficients are printed in ascending powers of S. 
ie The roots of the original and reduced equations 


are previded to enable the user to study the response cf the 


systems in the frequency domain. 
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See vetce Ore VariamlesS FOr print out in table form 
is available with up to eight variables maximum for any one 


cf three feossible runs. 


di. Choice of variables for graphical output with up 
to four curves may be plotted separately cr all cn one 


geaoh. 


EeecCEdyaredl  Olcpue Eesponse  t2 input excitations 
corsisting of a Step, Ramp, or Sinusoidal input are 
available. MThe graphs display the original output response 
compared to the lower order response and the error is 
pictted tc display the differences in response. Muitiple 


inputs may be used. 


Bee OUBRCUTINES 


IvwmeStrerOu tines are Utiiazed in ROUTHI. The rocts of 
the original and reduced equations are determined by 
Subroutine PRQD , wWaoaich was taken from [Ref 3], and the 
tables and graphs are determined by subroutine RSDUCT1 which 
moc Od| fication Of INTEG? ;, {Ref 3]. 


i> 





IV- FOURTH ORDER EXAMPLE 


Pe eo Lice PROCEDURE 


To illustrate the "Routh Approximation Method", A 


fourth order example was chosen for simplicity. 


The method, described in detail in [Ref 1], consists of 


three basic steps. The following transfer function is used 


to illustrate the procedure: 


nes t 


20 + 328 + 24% + 384 ¢ 


The first step is to compute what are termed Alpha and 


Beta coefficients from the Routh Table shown on page 19. 


—, 





Alpha, =0.625 pe al; 
Alpha,=1.684 6.3158 
Alpha,=3.008 a 
Alpha ,=6.3158 

Beta a; 





Beta, =0.03125 O. 25 
Beta,=0.0 re 
Beta.=0.03958 

Beta ,=0.0 
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Step two in the procedure is to obtain what are termed 
the Routh convergents, which are based on the following: 
Letting A, (s) and B, (s) denote the denominator and the 


numerator, respectively, of the eh Routh convergent, 1.ée., 


A, (s) = Alpha,sS + 1 

B, (s) = Beta, 

A, (Ss) = Alpha, Alpha,“ + Alpha.s ae 
B, (s) = Alpha,Beta,S + Beta, 


The general expression from [Ref 1] is the following: 


A, (Ss) = Alpha, (s)A,_, (s) + A, 9 (S) 
B,. (s) = Alpha, SB, _, (s) 7 By (Ss) + Beta, eee 7S aes 
with A_, (s)=0.0 B_, (s)=0.0 
Ay(s) =1.0 By (s) =0.0 


THe wROULhN Convergents for the fourth order example are 


the following: 
R, (s)= ORCel25 
Oo 2 ae 


R, (s)= 05056316 


1.05263S° + 1.68425 + 1.0 


es) = 0.158833S* — 0.0083 


A Gee aise Ge6aas + 1 


Ry (s)= 5? 


2087 7 32S° ot 249° + 8S + 1.0 


ny 





The third, and final step in the procedure is to apply 


what is termed a reciprocal transformation, defined by 


H(s) = =xR, (1/s) 


which 1s merely a reversal of the order of the polynomial 
coefficients. Thus, for the example given, the reduced order 
approximations are given by the following: 


H, (s) = OR03125 


S + 0.605 
H,(s) = 0.056316 

Seen INOS 6 3 
H,(s) = -0.0083S* + 0.158333 

s? + 3.6338* + 5.0665 + 3.166 
H,(s) = i; 


os cc gs? +2487 foes =) 20 


As expected, the fourth-order approximation is the same 
as the original equation. 

The poles of the approximants are illustrated in table 
IV.2. The poles of the approximants approach the dominant 
poles of the original higher-order equation as the order of 


the approximant is increased. 


18 















eae yl ON e eA 
ee Graeme 
u Uu 
Se Cee eae Gena ere 
eae oe eee See Seen 
uU u u u 
8 ae ae oe 
u u 
pu eu ane aoe 
as Lu’eudty—-2M=2M 
... Peudty— ema" Sm © yd tw—-°M=2M 
Mieae eS eydty— oP j= 2M Cmoeudty—2Po="m 
i. Sete Geen BOP encan enan BE en en ce ee 
Pp eps cP. TP 


2u/'n=" e908 
"u/'as* e109 


u 
CM/S eed 


uU 
$P5/'"5=leq0g 





2 /’m="eudtw 


Vu/ m=feudtw 


TysaP Cc 


D> Rc 


SP, /1P5_leudty 








ROUTH TABLE 


Table IV.1 


ng 








AS \R3 \ Re 


POLES OF ROUTH APPROXIMANTS 
Table IV.2 
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VeeeCOUPUTER SXANPLE 


hero eWee NTH CRDER SYSTEHW 


A seventh order system with the transfer function given 


by [Ref 2], as 


384x10/ 


G(s)=— Eee ri 3 ao 
s/443259°+626708°+36159008'+75114000$°+5539200008-... 


 .+1443200000S+384x10/ 


Opeln) CLaccorea form as 


384x10/ 


G(s)= 
($*42S+10) (S+10) (S+20) (S+80) (S+120) (S+200) 


was reduced to the low order approximants Dy 
utilizing ROUTH1 and the response to input excitations, 


consisting of a Step and Ramp are illustrated in figures 5.1 


wanouga 5.6. 


The rocts of the system and its low order approximants 


pee 2liusetrated in table ¥.1. 
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Poles and Zeros of Approximant 





-1.0 + 32.999 
R, |-10.89 + 32.30 {+ 316.58 
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ED wy 
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-2 

() , 
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POLES OF ROUTH APPROXIMANTS 
Table VV. 1 
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VeweecOrrantsON TO OTHER METACDS 


A. DOMINANT FOLE MeTHOD 


The Dentinant Pole Approximation method 1S paseG upon 
utilizing the poles closest to the imaginary axis. The 
eguation must be factored to obtain the characteristic 
rocts. Fer the given seventh-order system, the lower- order 


equations are given by this method as the following: 


10 

H (Ss) =—>————— 

S +2S5+10 

100 

SS ce 

S~+14S°~+30S+100 

2000 

al eee 

Seo oe oho +) UOS5+Z2000 

Graphical slots of the Dominant Pole reduced 


equations, in response to Step inputs, are illustrated in 
muagieecso,1 thorough 6.3 in comparison to the original 
equation and figures 6.4 through 6.6 illustrate tne 
ccmpariscn tc the Routh equations. The error between the 
systems is also plotted. Table VI.1 gives the analytical 


data for perficrmance measure comparisons. 
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Peet PR ALEVE O©PLTENTZATION METHOD 


The fcllowing equations, representing the seventh-order 
example, were taken from [Ref 2]. This method maxes use of 
an iterative finimization technique to locate the best fole 


and zero locations for the lower-order nodels. 


ac OS, 6 


Sa ss 
See eOUGet, «2056516 


piZecool 


S Z 


a 
Se eco oe toe oOloSoto2. 2861 


1470.1403 


~ 3 2 


H(s)=[>- 
ceepeo 2045 +208258425 +552.5241541470.1403 


Plots of these equations versus the Original 
equation and the Routh equations are illustrated in figures 
Segecntougn 6.12 in response to Step inputs. Table VI.1 


provides numerical data for performance measure compariscns. 
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See ReoENTATION OF DATA 


Comparisons were nade between the "Routh Approximation 
Method", and that of the Dominant Pole and Iterative 


Optimization methods previously described. 


iimeweasisesOm COMCarison consists of the following 


@riteria: 


1. Peak Overshoot = Ne” Va 


Maze Ss 


Poo hay liane -———-— > T,=time EOmmrce) tO reach 0.5 Y 


the very first time. 


S 


See. So ne-- —————— T =time Pome) co gomrcom O. 1 to 
0.9 of the final value. TH=1/5W. 


4. Settling Time----T_=time at which Y(t) =¥_ 
5. Graphical Representation in response to Step inputs. 


6. Average Error----J= IDE/t.,| ; t,=integration steps 


Table VI.1 illustrates the above comparisons. 
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R-3 aa 10.533 
: 7 wt 
; i : 
Eee 320:~C«~. 533°C 
| az : 
3 
p-3. | 1.346 [0.802 10.711 [5.866 1.15x10 “i 
| | | —— _ — 
pa | 1.083 10.400 Opera 15.866 0.062 | 
| as aes 
| | 
o-2 «jf 1.285 = 40.533 0.711 5.866 re | 
! E 
D-2 [1.350 )0-412 [0.684 |5.866 10.059 : 
TABLE VI.1 
Table A = Original Seventh Order 
symbols R = Routh Aporoximation 
D = Dominant Pole Approximation 
O = Optimum Minimization Method 
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Vitec OnCLUSTON 


The grafhicai data and tabulated analysis of this 
thesis irdicate that the Routh Approximaticn Methced is a 
valuakle formulation technigue which produces very 
satisfactcry results in acquiring approximations to 


higher-order systems with minimum cost. 


In ccmpériscn to the other methods cf analysis, the 
Routh ArEereximation Method offers a quick and easy 
analytical approach to obtaining lcw order models. Unlike 
the Fad'te Arpreroximation, the Routh method ensures statility 
of the lcwer-crder models, if the original higher-crder 
system is statle. The original system need net be factored, 


as in the Decninant Pcle Approximation method. 


The ccaputer program, ROUTH1, utilized to acquire the 
recuced crder equations, the roots cf the eguations and _ the 
graphical ficts and numerical tables, takes considerably 
less time than the minimization technique utilized in 
reference 2. However, the minimization téchnigue operates 
efficiently witheut prior knowledge of the system's transfer 
mane wl On. In the Routh method the availability cf the 


transfer function is a necessary requirement. 


In ccmparing the low order egquaticns of the various 
methcds discussed, the Routh method has proven to be a valid 
and efficient solution to the problem of cbtaining gocd low 


crder approxizants to complex higher-order systems. 
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